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ABSTRACT 

Observed clusters of galaxies essentially come in two flavors: non cool core clusters 
characterized by an isothermal temperature profile and a central entropy floor, and 
cool-core clusters where temperature and entropy in the central region are increasing 
with radius. Using cosmological rcsimulations of a galaxy cluster, we study the evolu- 
tion of its intraclustcr medium (ICM) gas properties, and through them we assess the 
effect of different (sub-grid) modelling of the physical processes at play, namely gas 
cooling, star formation, feedback from supernovae and active galactic nuclei (AGN). 
More specifically we show that AGN feedback plays a major role in the pre-heating 
of the proto-clustcr as it prevents a high concentration of mass from collecting in the 
center of the future galaxy cluster at early times. However, AGN activity during the 
cluster's later evolution is also required to regulate the mass flow into its core and 
prevent runaway star formation in the central galaxy. Whereas the energy deposited 
by supernovae alone is insufficient to prevent an overcooling catastrophe, supernovae 
are responsible for spreading a large amount of metals at high redshift, enhancing the 
cooling efficiency of the ICM gas. As the AGN energy release depends on the accre- 
tion rate of gas onto its central black hole engine, the AGN responds to this supernova 
enhanced gas accretion by injecting more energy into the surrounding gas, and as a 
result increases the amount of early pre-heating. We demonstrate that the interaction 
between an AGN jet and the IGM gas that regulates the growth of the AGN's BH, can 
naturally produce cool core clusters if we neglect metals. However, as soon as metals 
are allowed to contribute to the radiative cooling, only the non cool core solution is 
produced. 

Key words: galaxies: clusters: general - galaxies: active - galaxies: jets - methods: 
numerical 



1 INTRODUCTION 

Galaxy clusters are the largest virialised structures in the 
Universe directly observed using radio relics, the Sunyaev 
Zel'dovich effect, gravitational lensing, and X-ray measure- 
ments. These observations provide support to a cold dark 
matter like bottom-up formation scenario, by showing that 
a large variety of clusters, spanning an increasingly wider 
mass range, are present at different epochs, from smooth re- 
laxed objects to disturbed and merging ones. However, such 
support comes at a cost: as one goes back further in time, the 
build-up of clusters becomes a progressively more complex 
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process that involves accretion of cold material from the 
diffuse inter-galactic medium (IGM) and filaments, shock 
heating in the proto-cluster potential well along with multi- 
ple mergers potentially stripping a significant amount of gas 
and stars from galaxies. 

Although the assembly histories of different clusters 
are generally extremely different. X-ray measurements 
of unperturbed galaxy clusters feature gas with similar 
properties, especially entropjQ profiles, which are strikingly 
self-similar down to the central 100 (or even 10 in some 

^ We talce the usual definition of the entropy in astrophysics K = 
T/nl^^, where T is the gas temperature and is the electron 
number density assuming a mean molecular weight for electrons of 
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cases ) kiloparsecs jLloYd-Davies et al. 2000l: Piffaretti et al.l 
20051: iDonahue et al.l l2006l: iMorandi fc Ettoril l20o'7f 



Cavagnolo et all l2009l : ISanderson et all l2009l :l Pratt et alJ 
2010^ This is interpreted as the consequence of the for- 



mation of a large-scale gravitational-shock which heats 
the cold infalling gas from the IGM ( 10'' K) into an 
extremely hot and turbulent intra-cluste r medium (ICM) 
with temperatures around 10^-10^ K (|Tozzi fc NormanI 
l200lh . 

However in the cores (< 100 kpc) of these clusters there 
seem to exist a dichotomy in gas properties which has led ob- 
servers to split cluster samples into 'cool core' and 'non-cool 
core' clusters. Arguably the most relevant of these gas prop- 
erties is the entropy, as non-cool core clusters exhibit a well 
defined entropy floor in their centre, whereas in cool-core 
clusters the entropy profile decreases wit h decreasing ra- 
dius on all scales (,Morandi fc Ett ori 200j; ISanderson et al.l 
I2OO9I ). This lack of entropy floor in the central regions of 
some clusters is believed to be the signature of a cooling 
core because of the decreasing central temperature proflle 
associated to it fe anderson et al.ll2006l : IVikhlinin et al.ll2006l : 
IPratt et al., 2007 ). 

On the other hand, one of the long-standing problems 
of cosmological hydrodynamics simulations of cluster for- 
mation is the so-called 'cooling catastrophe' which takes 
place when gas is allowed to radiatively cool and 'standard' 
sub-grid physics, such as star forma tion and feedback from 
supernovae (SN), is implemented dSuginohara fc Ostrikei 



1998; iLewis et al 
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Ettori et al. 
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I2OOOI; 'Pearce et al 



2004; Kr avtsov et al 



200GI; 



20051 : 



Borgani et al 
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2OO7I ). Runaway cooling occurs, with se- 



vere consequences on both ICM and galaxy properties: the 
formation of an enormous cooling flow of gas (~ 1000 
M0 / yr) into the core triggers an anomalously high emission 
of X-rays, the presence of a persistent galactic disc of gas, 
and tremendous episodes of star formation. Central galaxies 
rapidly become too massive and blue. Thus, getting realistic 
thermodynamical properties of the ICM by numerical means 
constitutes a major challenge. 

Several processes have been invoked to solve this 
theoretical puzzle: pre-heating fr om early feedback pro- 
cesses in high-redshift galaxies ([Babul et al.l l2002h . Ac- 
tive Galactic Nuclei (AGN) fee d back from qu asars or 



radio jets jBinnev fc Taboij Il995l : iRephaeU fc sii3 Il995h . 
conduction of thermal energy fro m the outer shock- 
heated regions carried by electrons (|Voigt fc FabianI |2004| : 
iRasera fc Chandran 20081). and gas slosh ing from minor 
and major mergers "ijFabian fc Paine jll99ll ). Whilst the ver- 
dict is still out for early pre-heating and thermal conduc- 
tion, the ability of AGN feedback to stem cooling flows 
has already be en demonstrated in several hydrodynam- 
ics simulation s feiiacki fc Springelll2006': 'Siiacki et al.ll2007l: 
Khalatvan et a l. 2008; Pucliwcin ct al. 2008; Fabian et al.l 



2OIOI: iDubois et all |201GI : iMcCarthv et al.l I201ol . I2OIII : 
Teyssier et al.l 20111 ). It seems therefore natural to consider 



such a feedback mechanism as key to account for the self- 
regulation of the cold baryonic content of massive galaxies. 
However it remains to be demonstrated whether or not it 



fie = 1.167 in the hot cluster gas. The logarithm of this entropy is 
the standard definition of entropy in thermodynamics S = logK. 



can also explain the wide variety of observed ICM proper- 
ties. Especially, it is not obvious how such a physical process 
can reproduce the cool/non-cool core dichotomy observed in 
the population of gal axy clusters. 

In a recent work ([Dubois et al.|[201ol) . we demonstrated 
that jet-mechanical feedback from a central AGN is able to 
suppress the cooling catastrophe in a re-simulated cosmolog- 
ical galaxy cluster. Using a different numerical approach for 
gas dynamics and AGN feedback, [McCarthy et al.l ([201 1[ ) 
recently showed that the excess of entropy in the cores of a 
few galaxy groups is generated by the selective ejection of 
low-entropy material at high redshift z = 2-4 caused by both 
SN and AGN feedback. However much work is still needed 
to assess (i) whether the same mechanism works for more 
massive structures such as galaxy clusters, where mass loss 
from large outflows can be more easily prevented and (ii) to 
which extent such a mechanism depends on the numerical 
method and the specific subgrid model implemention of the 
AGN feedback. 

Therefore the aim of the present paper is to explore 
how various implementations of AGN feedback, which dif- 
fer either by the nature of the energy injected (kinetic or 
thermal) or by the epoch at which this energy injection oc- 
curs, impact the ICM gas properties. More specifically, we 
focus on the time evolution of thermodynamical quantities 
and mass distributions in our simulations, comparing and 
contrasting our results with observations whenever possible. 
The paper is organised as follows. In section [5J we describe 
the basic numerical ingredients of the simulations which are 
analysed. In section [3l we present our results and compare 
the evolution of entropy, density and temperature profiles of 
the simulated galaxy clusters. Finally in section|4l we discuss 
the consequences of these numerical experiments. 



2 NUMERICS 

The simulations presented in this paper use variations 
around the su b-grid physics and th e exact same initial con- 
ditions as in [Dubois et all ([2OI0I ). For completeness' sake 
we briefly recall in this section the basics of the numerical 
models that we employ. 



2.1 Physics of galaxy formation 

Gas radiates energy through atomic collisions a ssuming a 
H/He primordial composition ([Sutherland fc Etop ita 199 J). 
As a result it can collapse into dark matter (DM) poten- 
tial wells to form galaxies ([Silk[[ 19771 ). To model reioniza- 
tion, a homogeneous UV background h eating is imposed 
from z = 8.5 using the prescription of iHaardt fc Madaul 
([l996h . Star formation occurs in high gas density regions 
p > po {po = O.lH.cm"^). When the density threshold is 
surpassed, a random Poisson process spawns star cluster 
particles according to a Schmidt law p, — ep/ts, where 
is the gas free-fall time and e is the star formati on efficiency, 
taken to be e = 0.02 ([Krumholz fc Tan|[2007l ) i n order to 
repro duce the observational s urface density laws ( KennicuttI 
19981). The reader can consult iRasera fc Tevssieij ([20061 ) and 
Dubois fc Tevssi"eil ([20081 ) for more information on the star 
formation method. 
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Fe edback from SN ty p e II is included in a similar fash- 
ion as iDubois fc Tevssieil (|2008l ) : we model SN explosions 
by modifying the gas mass, momentum, and energy of sur- 
rounding cells following a Sedov blast wave solution. We 
adopt this approach because since we do not resolve SN rem- 
nants, thermal energy input in local cells would be quickly 
radiated away by the efficient gas cooli ng in the inter-stellar 
medium (ISM. lNavarro fc Whit3ll993l ). Our sub-grid model 
yields a physically motivated description of SN bubble ex- 
pansi on, and produces large-sc ale galactic winds in low-mass 
halos (|Dubois fc Tevssieij|2008l '). As we adopt a Salpeter IMF 
for star formation, we release an average of 10^^ ergs per 150 
Mq of stars formed in SN explosions. We return all the gas 
expelled by the SN explosion after 10 Myr, 10 % of which is 
assumed to be metals which we advect as a passive scalar. 
These metals also contribute to the cooling function of the 
gas assuming a Solar abundance ratio for the different ele- 
ments. 

We use a modified equation of state (EoS) at high gas 
density p > po to prevent numerical instabil ities from arti- 
ficial ly growing on the smallest grid scales (|Truelove et al.l 
More specifically, the minimum temperature in dense 
regions is set to T^in = To(p/po)"'\ with To = lO^'A', 
and n = 4/3 which leads to a constant Jeans mass A/j = 
1.3 10® Mq . Such a value of the polytropic index n is a rough 
approximation of the complex functional form of the EoS ob- 
tained by a nalytical modeling of the mu ltiphase structure of 
the ISM in lSpringel fc Hernguistl (|2003l ). 

We introduced a n ovel n umerical scheme for AGN feed- 
back in [Dubois et all (|2010l ) based on bipolar kinetic out- 
flows (or jets). We assume that a unique black hole (BH) en- 
gine located at the center of each galaxy pumps a fraction of 
the energy it gains by accr eting gas into such a jet. BHs, are 
modeled as sink particles (|Bate et al.lll995l : iKrumholz et al.l 
I2OO4I '). and can grow both by mergers with other BHs, and 
by accreting some of the surrounding gas according to the 
Bondi accretion rate 



Mi 



(1) 



where 01 = [p/ poY is a dimensionless boost factor that ac- 
counts for the unresolv ed small scale structure of the ISM 
ijBooth fc Schavell2009[ l. G is the gravitational constant, 
Mbh is the black hole mass, p is the average gas density, 
Cs is the average sound speed, and Ur is the average velocity 
of the black hole with respective to the surrounding gas. The 
accretion rate onto a BH is Eddington-limited 

47rGA/BH?Tlp 



A/] 



(2) 



where ax is the Thompson cross-section, c is the speed 
of light, mp is the proton mass, and is the ra- 
diative efficiency, assume d to be equal to 0.1 for the 
IShakura fc SunvaevI (| 19731 ) accretion onto a Schwarzschild 
BH. 

The total AGN luminosity is simply proportional to the 
rest-mass accreted energy 



-EaGN = e/ErA/BHC 



(3) 



where e/ is set to 1 to recover MeH- Aft, M-Bn-cy* relations 
consistent with obse rvational findings (|Tremaine et al.ll2002l : 
iHaring fc Rixl2004^ using AGN jets (Dubois et al., in prep.). 
The gas properties are locally modified around the BH to 



account for the presence of a jet with 10"* km/s vel ocity 
(as in the original prescription of lOmma et al.l (|2004h . see 
[Sb ois et al. I I2OIOI for more details on the scheme) within a 
cylinder of radius Aa; and height 2/S.x. 

We have al so introduced a the rmal implementation of 
AGN feedback ssier et al.ll201ll ) based on earlier works 
performed with smoothed-part i cle hydrodynamics (S PH) 
simulations (jSiiacki et al.ll2007l : [Booth fc Sch"av3l2009l ). In 
this alternative subgrid modelling, we release the AGN en- 
ergy in a thermal form as soon as it reaches a level high 
enough to reheat the surrounding gas to temperatures above 
10^ K. The radius of the injection bubble centred on the BH 
is chosen to be Aa;. As the energy is more efficiently coupled 
to the gas, we have to assume a different efficiency e/. As 
in the jet case, its value of 0.15 to reproduce the AfBH-Af*, 
A/bh-o"* relations. All other parameters of the model are 
kept identical to the AGN jet feedback. 

2.2 Initial conditions and simulation runs 

The simulations are ru n with the Ada ptive Mesh Refinement 
(AMR) code ramses (|Tevssieij [20021 '). The Euler equations 
are solved using a second-order unsplit Godunov scheme. 
More specifically, fluxes at cell interfaces are computed using 
an approximate Riemann solver where variables interpolated 
from cell-centered values are reconstructed with a first-order 
MinMod Total Variation Diminishing scheme. The Poisson 
equation is solved using a particle-mesh method where col- 
lisonless particles (DM, stars and sink particles) are added 
to the AMR grid with a Cloud-in-Cell algorithm. 

We assume a fiat ACDM cosmology with total mat- 
ter density Q.m ~ 0.3, baryon density fib ~ 0.045, dark 
energy density SIa = 0.7, rms mass fluctuation amplitude 
in spheres of 8/i~^Mpc ag = 0.90 and Hubble constant 
Ho = 70km.s~^.Mpc~^ that corresponds to the Wilkinson 
Microwave Anisot ropies Probe (WMA P) 1 year best-flt cos- 
mological model (|Spergel et al.l [20031 ). The simulations are 
performed using a resimulation (zoom) technique: the coarse 
region is a 128^ grid with Mom = 2.9 x 10^° M© DM resolu- 
tion in a SOh^^Mpc simulation box. This region contains a 
smaller 256^ equivalent grid in a sphere of radius 20 h~'^Mpc 
with AfDM = 3.6 X l(f Mq DM resolution, which in turn en- 
closes the flnal high resolution sphere with radius 6 h~^Mpc, 
a 512^ equivalent grid and A^dm = 4.5 x 10** Mq DM reso- 
lution. 

Cells in the highest resolution region may be adaptively 
refined up to ^max ~ 16 levels of refinement, maintaining a 
constant physical resolution of 1.2h~^kpc at all times. Re- 
finement follows a quasi-Lagrangian criterion: if more than 
8 DM particles lie in a cell, or if the baryon mass exceeds 8 
times the initial DM mass resolution, the cell is refined. A 
Jeans length criterion is also added to ensur e the numerical 
stability of the scheme on all levels £ < fmax (|Truelove et al.[ 
[l997h where 5p = p/p > 10^: we force the cells fulfilling 
these conditions to sample the local Jeans length with more 
than 4 cells. 

With this particular set of initial conditions we fol- 
low the formation of a galaxy cluster with mass Afsoo ~ 
2.4 10^* Mq, and radius rsoo = 940 kpc at 2 = 0, where the 
500 subscript stands for the mass and radius at 500 times 
the average density of the Universe (see AGNJETrun simu- 
lation data in table [1]). This cluster undergoes a 1:1 major 
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Figure 1. Time sequence of the galaxy cluster resimulation ZAGNJETrun at z = 4.03 (upper left panel), z = 3 (upper middle panel), 
z = 2 (upper right panel), z = 1 (bottom left panel), z = 0.6 (bottom middle panel), and 2 = (bottom right panel). This simulation 
includes gas cooling, UV background, star formation, SN feedback, metal enrichment, and jet AGN feedback. The gas density is shown 
in magenta, gas temperature in cyan, and gas metallicity in yellow. The image is a eight times zoom into the simulation box. 



merg er from z = 1 down to z = 0.6 (see fig. 2 of lDubois et al.l 
l20ld ). naturally separating early time evolution where tiie 
pre-lieating in the proto-cluster structures takes place, from 
late-time evolution during which the gas of the cluster re- 
laxes. 

In [Dubois et al.l (|2010l ) , we performed two simulations 
to explore the role of AGN feedback on the regulation of 
the cooling catastrophe: one simulation with AGN feed- 
back, cooling and star formation (AGNJETrun), and an- 
other identical to the first where AGN feedback was switched 
off (NOAGNrun). In this paper, we add several new simu- 
lations to this set to understand the impact of the AGN 
sub-grid modelling and the influence of early pre-heating on 
our results. Namely, we run an adiabatic simulation that 
ignores cooling and star formation ( ADIArun), a simula - 
tion with thermal AGN feedback as in lTevssier et al.l (|201ll ) 
(AGNHEATrun), and a simulation that involves restarting 
the AGNJETrun simulation at z = 0.58 (after the major 
merger has taken place) and thereafter deactivating the gas 
accretion onto the BHs and the related AGN feedback (AG- 
NOFFrun). The first new simulation (ADIArun) allow us to 
study the thermal properties of this cluster under the sim- 
plest physical conditions (gravity and gas dynamics only). 
The second new simulation (AGNHEATrun) allows us to 
compare the impact of the type of AGN feedback on the evo- 
lution of the cluster. Finally with the third new simulation 



(AGNOFFrun) we assess whether it is the AGN feedback at 
early or late times that shapes the thermodynamics of the 
cluster at z = and prevents a cooling catastrophe. 

On top of that, we also add a new set of simulations 
including metal cooling where the metals are produced self- 
consistently in SN explosions and advected with the fiow 
(ZNOAGNrun, ZAGNJETrun, and ZAGNHEATrun) . These 
enable us to quantify the importance of metal cooling on the 
ICM gas dynamics. 

The properties of the whole set of simulations are sum- 
marized in table [1] 

Fig. H] shows a time sequence of the ZAGNJETrun sim- 
ulation inside the zoom region (10 h~^.Mpc comoving on 
a side). The early evolution of the progenitors of the clus- 
ter at high redshift z > 1 (upper panels) shows the fila- 
mentary structure of the early Universe (gas density in ma- 
genta), and the deposit of metals (yellow color) in the IGM 
through SN feedback and AGN feedback, which correspond 
to a strong pre-heating phase. The merger of the two clus- 
ters between z = 1 (pre-merger phase, bottom left panel) 
and z = 0.6 (post-merger phase, bottom middle panels) pro- 
duces a strong shock (temperature coded in cyan) during the 
encounter and shock waves develop at larger distance during 
the post-merger. The z = final output of the simulation 
shows the relaxed cluster (bottom right panel). 
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Figure 2. Cuts of the gas entropy through the cluster core at z = for the NOAGNrun (upper left), the AGNOFFrun (upper right), 
the AGNJETrun (bottom left), and the AGNHEATrun (bottom right). Color bar units arc in log(keV.cm^). The picture size is 893 kpc. 



3 RESULTS 

We follow the time evolution of the entropy profiles, as well 
as other thermodynamical properties of the ICM, for our 
simulations without metal cooling (figure U]) and with metal 
cooling (figure [9]) in order to understand what drives the 
z — entropy profiles displayed in fig. [T] Each of the follow- 
ing sections systematically evaluates the consequences of in- 
cluding increasingly more physics into the simulations. First, 
in section 13.11 we map out the gas structure of the galaxy 
cluster at z — 0, and discuss how we separate the ISM of 
the central galaxy from the ICM gas. Section 13.21 presents 
results from the simplest case of cluster evolution under the 
influence of pure gravity and hydrodynamics (ADIArun). 
Section 13.31 considers the addition of atomic cooling, star 
formation, and UV reionization (NOAGNrun). Section [3.41 



studies the cluster evolution with the additional complexity 
of AGN feedback, comparing effects of two different subgrid 
models: AGN jet feedback (AGNJETrun) and AGN feed- 
back modeled as a thermal intput (AGNHEATrun). Sec- 
tion [3]5] assesses the role of pre-heating (AGNOFFrun), and 
finally section [331 examines how our results change with the 
inclusion of metal enrichment from supernovae (ZNOAGN- 
run, ZAGNJETrun, Z AGNHEATrun, ZAGNOFFrun). 



3.1 Azimuthal structure and radial profiles 

Fig. [2] shows a two-dimensional entropy cut through the cen- 
tre of the cluster at z = for the four simulations which 
include primordial gas cooling (NOAGNrun, AGNOFFrun, 
AGNJETrun, AGNHEATrun). It appears that away from 
the galaxy disk interface, the entropy distribution is smooth 
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Figure 3. Cuts of the gas temperature through the cluster core at z = for the NOAGNrun (upper left), the AGNOFFrun (upper 
right), the AGNJETrun (bottom left), and the AGNHEATrun (bottom right). Color bar units arc in kcV. The picture size is 893 kpc 



but non-uniform. Lower entropy gas is preferentially lo- 
cated in the central regions. Due to radiative losses from 
atomic cooling, a large disc component appears in two cases 
(NOAGNrun and AGNOFFrun), signaling a cooling catas- 
trophe (massive disc and high star- formation). This cold 
ISM component is also visible in fig. [3] with gas tempera- 
tures below 1 keV. In the simulations where AGN feedback 
proceeds down to z = (AGNJETrun and AGNHEATrun), 
the galactic disc is much smaller and the spatial extent 
of low entropy material is greater than in the cases where 
catastrophic cooling proceeds unimpeded (NOAGNrun and 
AGNOFFrun). The AGN jet in the AGNJETrun simulation 
manifests itself close to the disc by launching a bipolar out- 
flow of high entropy and high temperature ga s which sends 
soun d waves out to large radii in the ICM (see [Dubois et all 
I2OIOI ). 



Interestingly, the ICM in the NOAGNrun simulation 
has the highest level of entropy within the core, even though 
this is the simulation that does not include any feedback 
process at anytime. Compared to the other simulations, the 
central ICM also has a much larger temperature (fig. |3]). 
The AGNOFFrun simulation shows similar characteristics: 
the ICM develops a high-entropy, high-temperature phase 
on both sides of the disc. 

Fig. 31 and subsequent figures including 1-D radial 
gas profiles are obtained using volume-weighted angular- 
averaged profiles centered on the most massive galaxy at 
each time step of the simulation. The galaxies are identi fied 
using a halo finder (HaloMaker from iTweed et al.ll200^ us- 
ing the Most-Massive Sub-halo Method) on star particles or 
on DM particles when no stars are present. 

To fairly compare gas properties from simulations to 
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Table 1. Simulations performed with different gas physics and sub-grid models. Column 1 gives the name of the simulation, column 
2 indicates the form of cooling (none, H/He primordial composition, or metal), column 3 specifies whether star formation is activated, 
column 4 if supernovae are included, column 5 the form of AGN feedback, columns 6 and 7 list rsoo and M500 at z = respectively, 
columns 8 and 9 list r2ooo E'lid M2000 at 2 = respectively, column 10 lists the most massive BH mass Mbh in the central galaxy at 

2 = 0. 



Name 


Cooling 


SF 


SN 


AGN 


fsoo 


M500 


''2000 


A^2000 


Mbh 












(kpc) 


(Mq) 


(kpc) 


(M0) 


(Mq) 


ADIArun 


No 


No 


No 


No 


944 


2.3610^* 


486 


1.2910" 




NOAGNrun 


Primordial 


Yes 


No 


No 


956 


2.4610" 


514 


1.5210" 




AGNJETrun 


Primordial 


Yes 


No 


Jet CDubois et al. 2010) 


944 


2.3610" 


486 


1.2910" 


1.710^" 


AGNHEATrun 


Primordial 


Yes 


No 


Heat (Tevssier et al. 2011) 


949 


2.41 10" 


496 


1.3710" 


7.810"' 


AGNOFFrun 


Primordial 


Yes 


No 


Jet (2 > 0.58) & No (2 s: 0.58) 


946 


2.3810" 


496 


1.3710" 


1.1 lOi" 


ZNOAGNrun 


Metal 


Yes 


Yes 


No 


961 


2.5010" 


536 


1.7310" 




ZAGNJETrun 


Metal 


Yes 


Yes 


Jet (Dubois et al. 2010) 


939 


2.3310" 


474 


1.1910" 


3.3 10^" 


ZAGNHEATrun 


Metal 


Yes 


Yes 


Heat (Tevssier et al. 2011) 


936 


2.31 10" 


479 


1.2310" 


1.310" 


ZAGNOFFrun 


Metal 


Yes 


Yes 


Jet (2 > 0.58) &; No (2 s; 0.58) 


946 


2.3810" 


486 


1.2910" 


6.410^0 



observations (mostly from X-ray emission), one should use 
X-ray-weighted profiles instead of volume-weighted profiles. 
But the dense disc component in some simulations (AGN- 
JETrun and AGNOFFrun) pollutes the profiles even when 
cold cells (T < 1 keV) are removed. For this reason, we 
choose to only use volume-weighted profiles and remove gas 
cells with temperature smaller than 1 keV. We discuss the 
consequences of this choice in more detail in Appendix [XI 

3.2 The effect of gas dynamics 

The ADIArun simulation does not include any type of ra- 
diative cooling, star formation, or external heating process 
(i.e. no stellar or AGN feedback, and no UV background 
to model reionization), and thus serves as a reference for 
the more complex simulations. It tests the hydrodynamical 
evolution of the cluster under gravity. 

Entropy profiles in fig. [4] (first column, first row) show 
the well-known features of any adiabatic simulation per- 
formed with Godunov-type solvers: a self-similar power-law 
profile at large radii with a flat core component. The dotted 
line on the entropy plot corresponds to the r^'^ power-law 
inferred at large radii from both S PH an d AMR simulations 
of galax;y clusters by IVoit et al.1 (|2005l ). The high entropy 
material following a power-law distribution is believed to 
originate from a large-scale expanding shock at a few (~ 2- 
3) virial radii while the fiat core component develops from 
the lower entropy material, produced by small-scale turbu- 
lence and wea k shocks in the depths of the gravitational 
potential well (|Tozzi fc Normanll200ll ). Our results for the 
entropy profile in ADIArun are in very good agreement with 
the numerical study performed by I Mitchell et all (|2009l ). As 
pointed out by these authors, rather than producing a flat 
entropy core, adiabatic cluster simulations with standard 
SPH codes produce near power-law proflles down to the res- 
olution limit (smoothing-length) due to their inability to 
capture Kelvin-Helmoltz instabilities jAgertz et al.ll2007l ). 

The entropy profiles in the late time evolution spanning 
~ 6 Gyrs from z = 0.58 to 2: = for ADIArun vary only 
by ~ 20% due to minor mergers and gas turbulence. In- 
deed, the amount of entropy in the adiabatic simulation of 
this cluster is mainly determined by the major merger event 



lasting from z ~ 1 to 2 ~ 0.6 which doubles the entropy 
level of the pre-merger phase. Meanwhile, the gas tempera- 
ture increases as kinetic energy from the merger encounter is 
thermalised by shocks. The heating from this major merger 
could potentially prevent the formation of a massive cooling 
flow IZuHone et al.ll2010l ). However, even though radiative 
losses are not considered in the ADIArun, we can compute 
the theoretical cooling time based on primordial atomic rates 
(fig-Bl first column, fourth row). This cooling time increases 
from 400 Myr in the pre-merger phase at z = 1.2 to 700 Myr 
at the end of the merger z = 0.6. Thus the major merger 
would only slightly delay the occurrence of a cooling catas- 
trophe. The fact that the cooling time is still very small 
{t < 1 Gyr) after the major merger in comparison to the 6 
Gyr or so which remain to reach z = 0, means that several 
major mergers of this type would be required to balance the 
gas radiative losses. In the following sub-sections, we will see 
that this "sloshing" mechanism is not enough to prevent the 
cooling catastrophe. 

Like the entropy proflles, density proflles in the ADI- 
Arun (flg. [4] first column, second row) vary very little 
with time. Gas concentrates in a fiat core with low den- 
sity (~ 0.04 cm^''). In fact, the maximum density reached 
in the cluster core is lower than the gas density threshold 
for star formation (po ~ O.lH.cm""^). Thus to trigger and 
maintain star formation, the simulation needs to include a 
mechanism, like radiative cooling, which allows gas to effi- 
ciently concentrate in the gravitational potential well (see 
section |3.3|) . 

The temperature profiles observed in fig. [4] (first col- 
umn, third row) are consistent with gas in hydrostatic equi- 
librium and are therefore the direct conseque nce of the gas 
and total mass distribution in the cluster (see [Hansen et ahl 
I2OIII ). For gas in pure hydrostatic equilibrium, the pressure 
obeys: 

—f — = ~pG ^ '- , (4) 

ar r 

where p is the gas density, G the gravitational constant, 
r the spherical radius, and Mtot(< r) the total mass (gas, 
stars, DM and BHs) within r. Fig. [4] (first column, fifth 
row) shows that the gas pressure is 90% of the hydrostatic 



8 Y. Dubois et al. 



ADIArun NOAGNrun AGNJETrun AGNHEATrun AGNOFFrun 




10 100 1 10 100 1 10 100 1 10 100 1 10 100 

r (kpc) r (kpc) r (kpc) r (kpc) r (kpc) 



Figure 4. A comparison of tlie time evolution post major merger (z = 0.58) of the ICM volume-weighted, angular averaged entropy (first 
row), density (second row), temperature (third row), gas cooling time (fourth row), ratio of gas pressure to hydrostatic pressure (fifth 
row) and radial gas velocity dispersion (sixth row) profiles for the ADIArun (first column), NOAGNrun (second column), AGNJETrun 
(third column), AGNHEATrun (fourth column), AGNOFFrun (fifth column). Colors correspond to profiles at different redshifts listed 
on the entropy plots (first row). The dotted line in the entropy profiles corresponds to the r^'^ power-law inferred at large radii from 
both SPIf and AMR simulations of galaxy clusters by Voit et al. (2005). We have also plotted the pre-merger profiles of the cluster for 
the ADIArun (2 = 1.2) (dashed lines in left column panels). 



pressure. Some additional support for the cluster gas may 
be provided by the global rotation of the cluster and it s 
internal turbulence (jPolag et all l2005l : [Nagai et all l2007t ). 
To assess contributions from the latter, we measure the gas 
velocity dispersion (fig. [H first column, sixth row) in the 
inner parts of the cluster. It is 100—200 km/s, which is 5—10 
% of the sound speed in the relaxed (post major merger) 
phase and reaches 15% of the sound speed after the merger 
(at z — 0.58), in good agreement with recent simulations 
from IVazza et all (|2010l ). These values are extremely close 
to what is needed to make up for the lack of internal energy 
support inferred from the pressure profiles. 



We do not plot the radial velocity profiles for any of the 
simulations, as pressure profiles describe the dynamical state 
of the gas. However we comment that for the ADIArun, the 
radial velocity is negligible as it represents less than 1% of 
the sound speed in the cluster core. This demonstrates that, 
in the absence of any radiative cooling process, the gas is 
in (almost) perfect equilibrium in the core and is unable to 
concentrate in a dense galactic disc. 
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Figure 5. Cumulative gas (upper left), star (upper middle), baryon (logarithmic in upper right and linear in lower left), dark matter 
(lower middle) and total (lower right) mass profiles of the cluster at z = for the ADIArun (black solid line), the NOAGNrun (black 
dotted line), the AGNJETrun (black dashed line), the AGNHEATrun (black dot-dashed line), and the AGNOFFrun (green dashed line). 
The cumulative baryon mass is plotted on a linear scale in order to better display what is happening at large radii. 



3.3 The effect of atomic gas cooling 

When atomic cooling and star formation are allowed 
(NOAGNrun), gas can cool down to very low temperatures 
(10" K) compared to typical ICM temperatures (T ~ 10* K). 
Cold gas condenses in the centres of gravitational potential 
wells to form galactic discs. In turn, the gravitational po- 
tential deepens, pulling more and more DM into the centre. 
This well-known mecha nism called 'adiab a tic co ntraction' 
is described in det ail in iBlumenthal et all (|l986l ) (see also 
iGnedin et al.|[20o3 ). Displaying integrated mass profiles for 
the gas, stars, baryons and DM, Fig. [5] illustrates this ef- 
fect: as more baryons concentrate in the core of the cluster, 
more DM particles are pulled in, forming a stronger DM 
core component. 

Fig. U (second column, first row) shows that at all times, 
the entropy in the core of the NOAGNrun is one order of 
magnitude higher than the entropy in the core of the ADI- 
Arun. Furthermore, the entropy in the NOAGNrun continu- 
ously increases with time in the cluster, both inside and out- 
side the core. This entropy rise has been interpreted as aris- 
ing from the removal of low entropy gas as it cools and turns 
into stars (note that the cooling time of the ICM is shorter 
than 6 Gyr, the time span between z = 0.58 and z = 0) and 
its replacement with high- entropy material coming from th e 
outer parts of the cluster (|Brvanll200(il : IVoit fc Brvanll200ll ). 



Because of the Eulerian nature of grid codes, we cannot di- 
rectly follow the history of gas elements in our simulations. 
However using the Lagra ngian SPH techn i ciue, o r trac er par- 
ticles within grid codes, iMcCarthv et al.l (|201ll ) and IVazzal 
l|201ll ) have validated the above picture by showing that 
when low-entropy material condenses into the central disc. 



high-entropy material flows into the central parts of the ha- 
los replacing the depleted low-entropy gas. 

Support for the above picture also comes from the gas 
density profiles (fig. ID second row of second column) and 
the cumulative mass profiles (fig. [5]). The gas density in 
the NOAGNrun strongly diminishes with time in the ICM 
(fig. |4l second column, second row) as cold gas gets incorpo- 
rated into the galaxy. The small feature observed between 
30-40 kpc in the gas density and gas entropy profiles is char- 
acteristic of the flaring tail of the disc, which is dense, warm 
and low entropy, and thus difficult to separate from the 'true' 
ICM. Indeed, in Fig. [S] we plot the gas mass profile (dotted 
line) for all the gas (ICM and gas in the galaxy) in the 
cluster. Since we include the galaxy's gas, the cumulative 
gas mass profile for the NOAGNrun (dotted line) shows a 
strong contribution from the centre of the cluster. The outer 
regions (r > 50 kpc) in the cumulative gas mass profile are 
depleted of gas corroborating the picture that gas from the 
outer regions moves inward, cools, gets incorporated into 
the galaxy component and forms stars. Note that the to- 
tal baryon mass at rsoo — 1 Mpc is hardly affected (on the 
level of 10 %) by atomic gas cooling. Since the cooling time 
in these regions is larger than a Hubble time, the accretion 
of material on large scales is mainly driven by cosmological 
accretion, and not by radiative processes. 

More evidence for the above picture comes from pres- 
sure profiles (fig 21 second column, fifth row) They reveal 
that the gas in the NOAGNrun is far from hydrostatic equi- 
librium within r < 50 kpc resulting in a massive cooling 
flow (wr > 1000 km/s) towards the centre of the cluster. 
Since low-entropy material cools faster than high entropy 
material (fig. [4] second column, fourth row) , it continuously 
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flows onto the galaxy and is converted into stars. This re- 
moves pressure support at small radii, causing high-entropy 
gas at large radii to smoothly flow towards the centre as it is 
sucked into the cluster core. This explains the difference in 
the cumulative baryon mass seen at rsoo (table [1] and first 
panel in second row of fig. O for the NOAGNrun (black 
dotted) compared to the ADIArun (black solid line). 

Thus, we conclude that even though gas loses internal 
energy by radiative losses, the net effect of cooling is to 
replace low-entropy gas depleted onto the galaxy by high- 
entropy gas from large radii, and to fill the ICM with gas at 
higher temperature (see fig. |4l third column, second row). 
The rise in total mass (dotted curve in Aftot plot of fig. [5)) 
of the cluster core drives a gas internal energy increase to 
counterbalance the extra gravity. Consequently, the massive 
reservoir of internal energy at large radii supplies the cluster 
core with what is missing to support its own collapse: gas 
at higher entropy and higher temperature. 

The presence of this strong cooling fiow is also char- 
acterized by a vigorous, turbulent radial velocity dispersion 
(fig- HI second column, sixth row), but this cannot explain 
the increase in temperature as the radial velocity disper- 
sion remains very low compared to the gas sound speed 
(~ 10 — 20%). In a nutshell, we see all the features of a 
cooling catastrophe in the NOAGNrun: gas density is irre- 
versibly depleted in the ICM, which leads to an unrealistic 
amount of gas flowing towards the central galaxy, and larger 
values of entropy and temperature in the cluster (see also 
iNagai et al]|2007i ). 
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Figure 6. Time evolution for the net radial mass flux as a func- 
tion of radius for the AGNJETrun (upper panel) and the AGN- 
HEATrun (lower panel). Colors correspond to fluxes measured at 
different redshifts as listed in their respective panels. 



3.4 The effect of AGN feedback 

As was shown in [Dubois et al] (|2010l '). jet-induced AGN 
feedback solves many problems associated with the cooling 
catastrophe: galaxies are less massive, they form less stars 
especially at late-times, and the cooling-flow is moderate and 
self-regulated. In an effort to understand how t his works and 
what leads to the entropy proflle presented in [Dubois et al.l 
(|2010l ). we re-examine the AGNJETrun simulation. 

Following the cluster major merger a.t z — 0.58, the 
entropy in AGNJETrun shows a plateau in the core of the 
cluster at 100 keV.cm^ (fig-Hl third column, first row), which 
slightly decreases at late times as a moderate cooling flow 
develops and gas flows into the central galaxy. Neverthe- 
less despite the cooling fiow, high-entropy material does not 
replace the low-entropy material. The AGN regulates the 
amount of cold, low entropy material that gets incorporated 
into the galaxy by ejecting some of it back into the ICM, 
preventing high-entropy gas at large radii from flowing in in 
the process. Thus even though the cluster develops a small 
cooling flow, the amount of low entropy gas which flows into 
the centre is limited. 

Contrary to the supersonic cooling fiow in the NOAG- 
Nrun, that in the AGNJETrun is quiescent and largely sub- 
sonic {vr <^ Cs) as can be seen from the radial velocity dis- 
persion profiles (fig. m third column, sixth row). According 
to the cooling times plotted in figure|4] (third column, fourth 
row), the cluster core should experience a cooling catastro- 
phe in less than 2-3 Gyrs (a few 100 Myr in the very central 
< 10 kpc region) if no source of feedback is active. It ac- 
tually takes twice the time for the cluster to replenish the 



cold gas component in the galaxy and trigger a strong AGN 
activity at 2 = 0. 

Understanding the mechanism of jet-induced AGN feed- 
back is complicated by the variety of consequences it can 
have depending on the jet strength. In general, jet-induced 
AGN feedback stirs the cluster gas. Because the gas is com- 
pressible, kinetic energy imparted to the ICM can eventu- 
ally turn into thermal energy. There are two ways of doing 
this: weak jets can gently stir the ICM and increase its tur- 
bulence, or strong jets can violently shock the surrounding 
gas and reheat it. While weak shocks increase the thermal 
energy indirectly by producing a high level of turbulence in 
the core that cascades into more weak shocks that further in- 
crease the entropy, strong shocks directly act on the entropy 
by increasing the temperature. Unfortunately, disentangling 
these two modes is a difficult task that is beyond the scope 
of this paper. 

However, it appears that the effect of these different 
AGN jet modes on the ICM is to produce various en- 
tropy profiles in the cluster at different times. This could 
explain the large variety of entropy profil es in real X-ray 
clust e rs independent of their mass range (|Sanderson et al] 
I2OO9I : IPratt et all |201gI ) as these would mainly dependent 
on the AGN activity. Indeed with the AGNJETrun simula- 
tion, we are able to produce a cluster that exhibits different 
entropy profiles at different times: from a flat entropy core 
at z — 0.58 to a power-law entropy proflle at z — 0.0. 

The gas density profiles (fig. |4j third column, second 
row) lend support to the picture that AGN jet feedback 
stems the gas flow to the cluster core. Above an interme- 
diate radius (100 kpc), the gas density is extremely close 
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to the gas density in tiie ADIArun simulation, whereas in 
the NOAGNrun, gas is strongly depleted (gas density is four 
times lower at 100 kpc at z = in the NOAGNrun compared 
to the AGNJETrun). These conditions at intermediate radii 
and above hold over large time-scales, except at z = where 
a small cooling flow develops. The core's gas density grows 
with time due to the cooling flow, but because of the AGN 
jet feedback all this gas cannot cool and form stars as effi- 
ciently as in the NOAGNrun. As a result the ICM density 
also grows in the core of the cluster. 

The temperature profile (fig. 21 third column, third 
row) in the AGNJETrun simulation tends to an isothermal 
state at z = 0. After gas is shocked by the major merger 
(z = 0.58), the temperature is larger in the central 10 kpc 
of t he cluster. AGN act ivity is negligible during this phase, 
(see [Dubois et al]|20ld ) so it cannot explain the high tem- 
perature. Thereafter, gas cools diminishing the amount of 
thermal energy in the centre, and creating a small cooling 
flow around z = 0. 

Interestingly, a thermal AGN feedback mode (AGN- 
HEATrun) rather than a kinetic mode (AGNJETrun) cre- 
ates a high temperature core component. We interpret this 
as a consequence of direct thermal energy injection from the 
central AGN source which translates into a larger injection 
of entropy at small radii. By not directly acting on the in- 
ternal energy, the kinetic (jet) mode yields lower entropy 
profiles in the centre. 

Like the jet AGN feedback, thermal AGN feedback pre- 
vents the cooling catastrophe. Indeed, the cumulative stellar 
mass profile in fig. [S] shows that the amount of stars formed 
at z = is very similar for the AGNHEATrun and the AGN- 
JETrun. Gas mass and total baryon mass distributions are 
also very similar at intermediate and large (r > 25 kpc) 
cluster radii, therefore leading to similar ICM gas proper- 
ties on these scales (see fig. The only difference in the 
mass distributions appear in the core of the cluster: the total 
mass distributions (baryons plus dark matter) in the case of 
thermal energy and kinetic energy inputs are almost super- 
imposed, but differences of a factor of a few are present in 
the separate components in the inner 20 kpc. 

In fig. [51 we have represented the mass fiux of gas at 
different radii and different times within spherical shells 
centred on the cluster for the AGNJETrun and the AGN- 
HEATrun respectively. Negative and positive fiuxes compen- 
sate to give a very faint net flux in the central 10 kpc of 
the cluster (flg. E]) for both the AGNJETrun and the AGN- 
HEATrun. However at intermediate distances (r > 10 kpc), 
larger negative or positive net fluxes are measured, as a re- 
sult of wave propagation from the central AGN source. The 
net flux in the cluster outskirts is largely negative as the 
AGN feedback does not signiflcantly impact the gas at such 
large distances. 

Fig. [7] shows that the simulation that includes AGN jet 
feedback (AGNJETrun) all the way down to z = reaches 
the lowest entropy, with values below 10 keV.cm^ a few kpc 
away from the centre. In contrast, the simulation that has 
no feedback at anytime (NOAGNrun) has the largest en- 
tropy floor in the core {K ~ 400 keV.cm^). In flg. [3 we also 
give the observational entropy proflles of 20 X-ray selected 
clusters with cool cores (light blue lines) or no n-cool cores 
(light orange lines) from I Sanderson et al. (|2009l 'l in the tem- 
perature range 1-10 keV. The non-cool core clusters show an 
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Figure 7. Entropy profiles of the cluster at z = for the ADI- 
Arun (black solid line), the NOAGNrun (black dotted line), the 
AGNJETrun (black dashed line), the AGNHEATrun (black dot- 
dashed line), and the AGNOF Frun (green dashed line ). Entropy 
profiles from real clusters from lSanderson et al. I boool ') with cool 
core clusters in light blue, and non-cool core cluster in light orange 
are overplotted on top of our simulation results. 



entropy floor in the centre of the cluster (light orange lines), 
whereas cool clusters possess entropy proflles declining to- 
wards the centre (light blue lines). Thus, the presence of jet 
AGN feedback allows the formation of a cool-core cluster 
whereas with the AGN heating mode the cluster is always 
maintained in a non-cool core conflguration. It is worth to 
note that without AGN feedback the entropy floor is kept 
very high in the core with an amplitude higher than those 
obse rved for the most ma ssive of the non-cool core clus- 
ters (|Sanderson et aLlbOOgl ). 

Finally, one might wonder how robust our proflles are to 
the choice of averaging method (spectroscopic-like/X-Ray or 
volume weighted) for the physical quantities. As discussed 
in detail in Appendix A, large differences between aver- 
ages appear in the core of the cluster (inner 50 kpc or so) 
when a massive disc of cold gas has developed in the cen- 
tral galaxy, fed by an unchecked cooling catastrophe. Oth- 
erwise, the agreement between all estimates is excellent, as 
demonstrated by the z > 0.40 proflles (black and dark blue 
curves in all panels) on Fig. lAll Therefore, at least all the 
inner proflles we measure in the (more realistic) runs where 
the cooling catastrophe is prevented (AGNJET and AGN- 
HEAT) are rehable. 



3.5 The effect of the pre-heating 

In the interest of understanding whether a strong initial pre- 
heating of the gas is alone able to prevent a cooling catas- 
trophe for several Gyrs, we run the following simple experi- 
ment: we restart the AGNJETrun simulation after the major 
merger, i.e. from z = 0.58, deactivating BH accretion and 
AGN feedback. We call this simulation the AGNOFFrun. In 
the literature, the term 'pre-heating' is generally employed 
in a generic sense to describ e a rise in the level of entropy 
of th e proto-cluster gas (e.g. lKaiseilll99il : lEvrard fc Henrvl 
which can be caused by various physical processes. 
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Figure 8. Time evolution as a function of radius of the net ra- 
dial mass flux for the AGNOFFrun. Colors correspond to fluxes 
measured at different redshifts as listed in the panel. 

In the current work, we refer to pre-heating as the rise in 
entropy which is self-consistently provided by the feedback 
processes modelled in the simulations before the final as- 
sembly of the cluster. This explains why the end of the pre- 
heating phase coincides with the beginning of the last major 
merger undergone by the cluster at 2 ~ 1. 

Between z = 0.58 and z = 0.42, 2 Gyrs of evolution 
elapse. Fig. 3] (fifth column, fourth row) shows that this co- 
incides with the cooling time for AGNOFFrun's ICM gas at 
z = 0.58 and r = 20 kpc. Therefore as expected, inside this 
radius, the entropy profile strongly declines during this time 
interval (fig. [4] fifth column, first row), signaling a cooling 
flow into the centre. More specifically the entropy in the clus- 
ter core (r < 10 kpc) decreases from its values a.t z — 0.58 
{K ~ 20keV.cm^) to very low values {K < lOkeV.cm^) at 
z = 0.42. Concomitantly the gas density in the centre of the 
cluster increases (fig.[ll fifth column, second row) as internal 
energy is radiated away by gas cooling, leading to a loss of 
pressure support (fig. [4] fifth column, fifth row). Neverthe- 
less by redshift z = 0.42, the cluster has not yet reached the 
cooling catastrophe stage; the cooling fiow is still moderate 
and subsonic. 

Gradually after z — 0.42, the cooling catastrophe devel- 
ops: gas flows supersonically into the gravitational potential 
well (see radial velocity dispersion profiles in fig. [l] fifth col- 
umn, sixth row), concentrates in the cluster centre (fig. [SI 
left panel in upper row) and forms a cold disc component 
(see fig. [3]). As the gas condenses in the centre, strong star 
formation occurs in the central galaxy (see cumulative stel- 
lar mass profiles in upper middle panel of fig. [S]) resulting in 
a more massive galaxy at 2 = in the AGNOFFrun than in 
the AGNJETrun. DM also concentrates more in the cluster 
centre through adiabatic contraction (fig. [5] middle panel 
in bottom row). Bizarrely enough, it is also clear from this 
figure, that the AGNOFFrun cluster features a DM halo 
less concentrated than in the AGNHEATrun. This, in turn, 
demonstrates that a different implementation of AGN feed- 
back can have a significant and counterintuitive impact on 
the 2 = DM concentration of clusters: the mode that re- 
duces the baryonic mass the most in the cluster core (AGN- 



HEATrun) at 2 = is also the mode which has the largest 
DM halo concentration. This is because the stellar compo- 
nent should be taken into account when calculating adia- 
batic contraction as it responds in the same coUisionless way 
as the dark matter. In other words, when one compares the 
total mass of dark matter plus stars (fig. [SI right panel in 
bottom row, at least in the center where the gas mass in neg- 
ligible), the AGNOFFrun becomes more concentrated than 
its AGNHEATrun counterpart, as naively expected. 

We can monitor the effect of preheating on the progres- 
sion of the cooling catastrophe in fig.[Hl It shows that as soon 
as the AGN feedback is turned off, the low-entropy mate- 
rial quickly fiows into the centre and becomes part of the 
central galaxy ISM. Gas located at larger distances is then 
subsequently accreted. In other words, higher entropy ma- 
terial replaces low-entropy gas thereby increasing the ICM 
entropy in the core of the cluster (fig. [4] fifth column, first 
row). In the meantime, fig. [4] (fifth column, first row) shows 
that the entropy continues to grow at larger radii (r > 50 
kpc) because of the loss of pressure support at increasingly 
larger radii as low entropy gas is removed from the clus- 
ter centre to form stars. Unlike the supersonic flow at small 
radii, at large radii (r > 100 kpc) the gas fiows quiescently 
into the core (see gas radial velocity dispersion in fig. [4] fifth 
column, sixth row). We also note that the entropy profile at 
2 = is larger in the core (r < 50 kpc) than at moderate 
distance (100 kpc) from the centre (see fig. [31 fifth column, 
first row). This cannot be the consequence of an adiabatic 
process, because the entropy profile would be monotonic. 
Therefore, there are two possible explanations for the be- 
havior of the 2 = entropy profile: (i) either it is the result 
of spurious entropy produced by the grid code which is not 
a strictly entropy-conserving scheme, or (ii) a strong shock 
developed somewhere in the cluster core causing the gas 
to endure a non-reversible process. Even though we cannot 
entirely rule out the former possibility, there is plenty of 
evidence in favor of the latter. As the central ICM rapidly 
flows onto a cold and dense galactic disc component, there is 
a strong discontinuity both in gas density and temperature 
at the disc interface. This is evidence of a strong shock. As 
a result, upwind of the shock the ICM is heated, and both 
temperature and entropy increase. Figs. [2] and [3] show this 
high-entropy, hot region perpendicular to the galactic disc. 



Finally, an interesting feature of the AGNOFFrun sim- 
ulation is that compared to the AGNJETrun which includes 
kinetic AGN feedback all the way down to 2 = 0, the gas 
properties at large radii are hardly modified by the absence 
of late-time AGN activity. This is because the cooling flow 
does not extend to radii greater than 100 kpc and the late 
AGN jet activity hardly perturbs the gas beyond these dis- 
tances. In light of our previous analysis, we conclude that 
pre-heating from AGN feedback at high redshift (2 > 1) can- 
not single-handedly prevent a cooling catastrophe, but that 
it plays an important role in (i) setting up the gas properties 
at large distances from the centre, and (ii) to some extent 
softens the consequences of the catastrophic cooling in the 
core (lower entropy, temperature of the gas and even mass 
of dark matter and baryons). 



Cluster entropy profiles 13 



ZAGNHEATrun 



ZAGNOFFrun 




' (kpc) 



- (kpc) 



' (kpc) 



■ (kpc) 



Figure 9. A comparison of the time evolution post major merger (z = 0.58) of the ICM volume-weighted, angular averaged entropy (first 
row) , density (second row) , temperature (third row) , gas cooling time (fourth row) , ratio of gas pressure to hydrostatic pressure (fifth row) 
and radial gas velocity dispersion (sixth row) profiles for the NOAGNrun (first column), ZNOAGNrun (second column), ZAGNJETrun 
(third column), ZAGNHEATrun (fourth column), and ZAGNOFFrun (fifth column). Colors correspond to profiles at different redshifts 
listed on the entropy plots (first row). The dotted line in the entropy profiles corresponds to the r^'^ power-law inferred at large radii 
from both SPH and AMR simulations of galaxy clusters by Voit et al. (2005). 



3.6 The effect of metal cooling 

Until now, this paper has discussed the results from the 
cluster resimulations without supernovae and the accom- 
panying stronger cooling due to metals released in the ex- 
plosions. In this section, we analyze versions of the cluster 
simulations with supernovae (ZNOAGNrun, ZAGNJETrun, 
ZAGNHEATrun, ZAGNOFFrun) and focus on the effect of 
metal cooling on the gas properties of the ICM. 

Like Fig. [4] for the no metal case. Fig. [9] shows the time 
evolution of different thermodynamical gas profiles for the 
ICM when SN and hence metal enrichment is included in 
the simulations. For ease of comparison, we display the re- 
sults from the simulation with primordial cooling (NOAG- 



Nrun) in the same panel. In the absence of AGN feedback 
(ZNOAGNrun), the galaxy cluster endures a strong cooling 
catastrophe as evidenced clearly by the gas pressure profiles 
(Fig. [9l second column, fifth row), even though energy from 
SN explosions is deposited into the surrounding gas. This is 
not a surprising result as it has been known for some time 
that the energy from SN feedback is not enough to suppress 
the c ooling catastrophe in massive structures (Nag ai ct al] 
I2OO7I ). Cosmological re-simulations of the formation and 
evolution of a proto-galactic structure i n a Milky- Way like 
halo performed at sub-parsec resolution (| Powell et al.ll201ll ) 
show that the mass-loading factor of large-scale, supernova 
driven outflows is low and therefore they cannot remove 
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Figure 10. Cumulative gas (upper left), star (upper middle), baryon (logarithmic in upper right and linear in lower left), dark matter 
(lower middle) and total (lower right) mass profiles of the cluster at 2 = for the ADIArun (black solid line), the NOAGNrun (black 
dotted line), the AGNJETrun (black dashed line), the AGNHEATrun (black dot-dashed fine), the ZNOAGNrun (blue dotted fine), the 
ZAGNJETrun (blue dashed line), and the ZAGNHEATrun (blue dot-dashed line). The cumulative baryon mass is plotted on a linear 
scale in order to better display what is happening at large radii. 



large amounts of gas from high-r edshift galaxies, as was als o 
shown in an idealized context bv lDubois fc Tevssieij (|2008h . 

Therefore, instead of alleviating the cooling catastro- 
phe, SNe actually produce a stronger one, with the total 
mass being more concentrated in the centre of the cluster 
than in the absence of SN feedback (NOAGNrun, fig. [TO)) . 
This is because the metals released by SNe reduce the cool- 
ing timescales at high redshift. Indeed, even though this ef- 
fect is very pronounced at lower temperatures (T< lO^K), 
and more moderate for the ICM (20 % variation at T = 3.5 
keV), the high redshift progenitors of the cluster are filled 
by material at much lower temperature than a few keV, 
and so metals lead to an increase of the gas cooling effi- 
ciency by about an order of magnitu de or more in these 
structures (jSutherland fc Dopitalll993l ). It results in larger 
cooling times at low redshift (see cooling timescale profiles 
in Fig. 1^ fourth row, first and second columns) because the 
cluster has suffered more from the cooling catastrophe and 
the gas that stays in the core of the ICM has a larger entropy. 

In other words, in the simulations with metals, a larger 
amount of gas collapses early on into galaxy and group size 
dark matter halos and forms stars. Halos thereby endure a 
more marked depletion of their gas content than in the com- 
parable simulation with no metals (NOAGNrun) , leading to 
lower gas density profiles for ZNOAGNrun (see Fig. |9] sec- 
ond row, first and second columns) . As explained in previous 
sections, as central regions get depleted, gas from the outer 
regions can flow into them because of the lack of pressure 
support. Because this loss of pressure support associated 
with the cooling catastrophe in the ZNOAGNrun not only 
takes place earlier but also at larger radii than in the NOAG- 



Nrun, gas flows into the cluster centre from larger distances 
in the ZNOAGNrun. Since the gas density is lower at large 
radius but the temperature of the infalling shocked IGM gas 
is the same as at smaller radii, the amount of entropy car- 
ried by gas coming in from larger radii is larger. As a result, 
higher entropy gas than that in the NOAGNrun replaces 
the low entropy material at small radii in the ZNOAGNrun, 
leading to entropy profiles that rise with time to a higher 
level in the run with metals (Fig. |9] first and second col- 
umn, first row). 

The simulations ZAGNJETrun and ZAGNHEATrun, 
which include AGN feedback with the kinetic and thermal 
modes respectively, show very similar gas properties (Fig. [Oj 
third and fourth columns), and mass distributions (Fig. llUp . 
Even the concentration of the cluster DM profile which was 
noticeably different in the AGNJETrun and AGNHEATrun 
simulations is now very similar. Thus, we will not consider 
hereafter the different effects of these two modes on the 
gas, which are already described for their equivalent sim- 
ulations with no SN (AGNJETrun and AGNHEATrun) in 
section l374l Instead we focus on how metals alter the effect of 
AGN feedback. First we compare the three simulations with 
metals, where one has no AGN feedback (ZNOAGNrun) and 
the others do (ZAGNJETrun, ZAGNHEATrun). As was the 
case without metals, for the simulations with metals (ZA- 
GNJETrun and ZAGNHEATrun) AGN feedback diminishes 
the amount of entropy in the ICM compared to the simu- 
lation without AGN (ZNOAGNrun). AGN feedback, again, 
suppresses the cooling catastrophe, even though more effi- 
cient metal cooling at early times makes this more challeng- 
ing than in the no metal case. We emphasize that, as in 
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the no metal case and despite a much larger amount of pre- 
heating from SN and AGNs at high redshift, only the runs 
where AGN feedback is present during the whole evolution 
of the cluster are able to prevent the cooling catastrophe. In 
the ZAGNOFFrun simulation (Fig. |9] last column) where 
AGN feedback is switched off after the last major merger, 
the central galaxy ends up accreting cold gas at a rate of 
« 500 Mq/ji by 2; = 0. All mass components (gas, stars, 
DM) at 2 = in ZAGNJETrun and ZAGNHEATrun are 
less concentrated in the centre of the cluster compared to 
the ZNOAGNrun (Fig.[10ll, and unlike in ZNOAGNrun, the 
gas supports itself against its own collapse in ZAGNJETrun 
and ZAGNHEATrun (Fig. [9l third and fourth columns, fifth 
row). 

Next we examine the differences between the AGN feed- 
back runs with and without metals. The presence of metals 
plays a crucial role in how the AGN impacts the gas. Com- 
paring entropy profiles at 2; = for simulations with AGN 
feedback and metals (ZAGNJETrun and ZAGNHEATrun), 
to the equivalent simulations without metals (AGNJETrun 
and AGNHEATrun) in Fig. [Til we see that, in the simu- 
lations with metals, entropy levels in the cluster core are 
larger by almost one order of magnitude. The cause of this 
effect can be traced back to the fact that metals increase the 
cooling efficiency of the gas in the early phases of structure 
formation. This leads to more fuel for AGN activity at early 
times which pushes more gas to larger distances from the 
centre of the halos. Note however, that this is only an inte- 
grated time effect. As a consequence, very early on, before 
the AGN has reached a mass when its feedback allows it to 
significantly impact the fate of the intra-halo gas, baryons 
are more concentrated in the ZAGNJETrun than in the AG- 
NJETrun. This is reflected in the DM density profile which 
is more concentrated in the ZAGNJETrun than in the AG- 
NJETrun (fig. [10] middle panel of bottom row), even though 
the amount of baryons in the ZAGNJETrun at 2 = is 
much smaller than in the AGNJETrun. 

Nevertheless, in general, because of these early time ef- 
fects, the integrated mass of baryons at 2 = at all radii 
(Fig. nop is smaller when metals and AGN are present (ZA- 
GNJETrun and ZAGNHEATrun) than for the equivalent 
AGN simulations without metals (AGNJETrun and AGN- 
HEATrun). The total stellar mass at 2 = in the cluster 
is also reduced by the combined effect of AGN activity and 
metal release by SNe. This underlines the fact that AGN 
feedback has a stronger effect on the cold baryon content, be- 
cause of the enhanced cooling rates which allow the BH cen- 
tral engine to grow faster at high redshift and thus quench 
star formation earlier on. This is also reflected by the fact 
that the central BH is more massive when metal enhanced 
cooling is enabled (see table [1} and implies larger amounts 
of energy deposited into the ICM. 

We note that the main results we obtain when metal 
cooling and AGN feedback are included in the simulations 
(i.e. increased amount of pre-heating at high redshift, re- 
duced stellar mass for the central galaxies, too high level 
of entropy in the cluster core, central temperature profiles 
which are too steep) are broa dly consistent with those ob- 
tained bv lFabian et al.l (|2010l ) for different clusters using a 
different numerical technique and a slightly different model 
of AGN thermal feedback. We also remark that we still 
have an adiabatic contraction of the DM halo, contrary to 
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Figure 11. Entropy profiles of the cluster at 2 = for the ADI- 
Arun (blaek solid line), the NOAGNrun (black dotted line), the 
AGNJETrun (black dashed line), the AGNHEATrun (black dot- 
dashed line), the ZNOAGNrun (blue dotted line), the ZAGN- 
JETrun (blue dashed line), and the ZAGNHEATrun (blue dot- 
dashed line). 



what [rtevssier et al.l (|201ll ) obtained in their Virgo-like clus- 
ter with thermal AGN feedback. In this paper, we inject 
the energy of the thermal feedback on the same scale as 
the jet mode, close to the resolution limit, to facilitate a 
comparison between the two ru ns. Note that t his is differ- 
ent from the scale adopted by iTevssier" et al.l lioil), who 
injected the thermal energy into large bubbles, and explains 
why we are not able to halt adiabatic contraction and they 
are. In short, injecting thermal energy on larger scales in- 
creases the effective efficiency of AGN feedback. A detailed 
study of parameter/resolution impact on our results is de- 
ferred to a companion paper (Dubois et al, in prep). 

At late times 2 < 0.58, there is no longer a cold gas com- 
ponent in the centre of the cluster of ZAGNJETrun and ZA- 
GNHEATrun (Fig. m third and fourth columns, third row) 
and the gas remains in hydrostatic equilibrium (Fig.|9] third 
and fourth columns, fifth row) because cooling times are 
comparable to the last major merger times (Fig.[9l third and 
fourth columns, fourth row), and therefore cold gas only trig- 
gers some faint AGN feedback activity. Note however that 
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such a reduced level of AGN feedback activity at low redshift 
is key to prevent a cooling catastrophe, as demonstrated 
by the gas profiles measured in the ZAGNOFFrun simula- 
tion (Fig. [9l fifth column) which lacks such a low redshift 
AGN feedback activity. The strong enhanced early phase of 
AGN activity in ZAGNJETrun and ZAGNHEATrun leads to 
lower entropy compared to the no metal simulations AGN- 
JETrun and AGNHEATrun at large distances from the cen- 
tre, because of the redistribution of the gas from the centre 
of the cluster to the outskirts (Fig. Illl bottom panel). How- 
ever, the temperature profiles for ZAGNJETrun and ZAGN- 
HEATrun seem rather unrealistic, when compared to obser- 
vations: they continuously increase with decreasing radius 
with extremely high values (iO-20 keV) in the central 10 
kpc (Fig. [21 third and fourth columnds, third row) . This is 
a direct consequence of the mass distribution in the cluster 
and the hydrostatic equilibrium. It is possible that a more 
complete AGN feedback prescription such as one that in- 
cludes photo-heating cou ld, at least partia lly, alleviate this 
problem. As suggested bv lCantalupol (|2010l ), soft X-ray pho- 
tons emitted by stars during strong starburst episodes or by 
AGN, can offset radiative losses by photo-heating the main 
ion coolants in proto-cluster and massive galaxy halos. If 
this kind of feedback takes place at high redshift, it could 
halt the cooling catastrophe without removing much mass 
from the centre. This would increase the central gas density, 
decrease the 'apparent' cooling time (as opposed to the effec- 
tive cooling time, which is increased by the photo-ionization 
effect), and allow for lower temperatures in the cores of the 
clusters. We defer a numerical investigation of this scenario 
to future studies. 

Finally, we must also point out that we measure a lower 
gas metallicity, Z ~ 0.15 Zq z = Q, plotted in fig. 1121 and 
a steeper s patial gradient that t hose observed {Z ~ 0.3- 
1 Zq from I Sanderson et al.l l2009t ) . This is even worse for 
the stellar metallicity, which i s observed to be lar ger than 
solar for giant ellipticals (e.g. ICasuso et al.l (|l996l )). Aside 
from uncertainties in the stellar yields and the stellar ini- 
tial mass function, there are (at least) two other reasons 
why we are not able to achieve more realistic metal abun- 
dances. The first one is the inability at a given DM mass 
resolution to resolve well enough the star formation of small 
satellites galaxies that are responsible fo r the enrichment of 
the I GM with large-scale galactic winds (|Dubois fc Tevssied 
l2008l 'l. The second is the lack of modeling of both stellar 
winds and Type la SN that return a large fraction of their 
stellar mass content back to the gas. We plan to address 
both issues in a near future but the comparison of the no 
metal and (low) metal runs that we have performed strongly 
suggest that a higher amount of metals in the gaseous phase 
will simply lead to more intense early AGN feedback. 



4 CONCLUSIONS 

We have performed high-resolution resimulations of a cos- 
mological galaxy cluster, using standard implementations 
for galactic physics (e.g. radiative cooling, star formation, 
supernova feedback, uniform UV background heating to 
model rei onization) and different prescriptions for AGN 
feedback (jPubois et al.ll2010l : pRvssier et al.ll201lh . 

Whereas observations indicate that non-cool cores are 
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Figure 12. Stellar metallicity profiles (first row) and gas metal- 
licity profiles (second row) for the ZNOAGNrun (left column) and 
the ZAGNJETrun (right column). 

high ly correlated with the level of disturbance of the ICM 
(e.g. [Sanderson et ahllioogl : iRossetti fc Molendill2010l ). sug- 
gesting that major mergers are responsible for the forma- 
tion of non-cool core objects, our set of simulations leads us 
to conclude that merger activity alone (i.e. not backed by 
AGN feedback) cannot sustain core profiles over timescales 
comparable to the Hubble time. It therefore seems difficult 
to argue that this mechanism can account for a significant 
fraction of the non-cool core cluster population. These simu- 
lations also show that pre-heating by SN and AGN feedback 
at high redshift alone, before the cluster has fully assembled 
is not sufficient to prevent the occurrence of a cooling catas- 
trophe. Instead, they lend support to the view that AGN 
must provide some amount of extra pressure support to the 
gas in the cluster core after the galaxy cluster has formed. 
This late feedback is also key to shape the thermodynamical 
properties of the central ICM. 

Indeed, we established that for a gas with pristine com- 
position, the natural interaction between an AGN jet and 
the ICM which regulates the growth of the AGN's BH, can 
produce entropy profiles with different shapes: sometimes 
the profile has an entropy floor in the core (i.e. the cluster 
will be classified as non-cool core), while at other times, if 
cooling is sufficiently rapid there will be a power-law en- 
tropy profile extending from large radii down to the cluster 
centre (i.e. the cluster will belong to the cool core cluster 
category). Hence, AGN jet feedback appears to be a good 
candidate to explaii i the variety of entropy profiles observed 
in galaxy clusters ([Sanderson et al]|2009l ). provided metal 
cooling is neglected. 

Without significant AGN feedback at intermediate and 
low redshift, clusters rapidly transfer their low-entropy ICM 
material to a galactic disc and replace it with the high- 
entropy gas stocked in their outskirts. As a result, a high- 
entropy fioor and large amounts of baryons are systematic 
features of cluster cores. Moreover, as the cooling catas- 
trophe ensues, gigantic starbursts erupt within the mas- 
sive galactic discs anchored at the bottom of the grav- 
itational potential well of these clusters. Taken at face 
va lue, these results seem to contradict the recent findings 
of [McCarthv et akl (|201l[ )'s simulations of galaxy groups. 
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In these lighter structures, they conclude that, on aver- 
age, turning off the AGN feedback at intermediate redshift 
(0.375 < z < 1.50) lowers the entropy of the gas in the core 
compared to simulations where AGN feedback is permitted 
to proceed down to z = 0. Assuming that the differences 
between our results cannot be attributed to numerics (SPH 
versus Eulerian grids) , we interpret them as an indication of 
a fundamental difference in the way feedback mechanisms 
can operate in galaxy groups and clusters. For the former, 
the effects of high redshift AGN feedback are still felt at 
z = 0, whereas for the latter they are wiped out by the 
stronger gravitational forces on time scales comparable to 
the assembly time and therefore need to be sustained at 
later epochs. 



Finally, we assessed the effect of metal cooling on the 
gas properties of the ICM. Metals ejected by SNe accelerate 
the gravitational collapse of baryons by sapping their inter- 
nal energy. As a result, in the absence of an AGN feedback 
mechanism, the cooling catastrophe is exacerbated. How- 
ever, when AGN feedback is turned on, AGN activity is 
fueled by a more rapid growth of the central BH engine, 
and thus has a more dramatic impact on the ICM and the 
baryon content of galaxies. Indeed as the gas concentrates in 
structures earlier, larger BH accreting at a higher rate can 
push the gas further away from the centers of less massive 
halos, resulting in less concentrated gas and stellar profiles in 
the final galaxy cluster at low redshift. Since the gas density 
decreases, and entropy scales as K (x rie^^^ , larger entropy 
floors are logically present in the core of metal rich galaxy 
clusters. These effects systematically lead to the formation 
of non-cool core clusters regardless of how AGN feedback 
was implemented in the simulations, which is in obvious 
contradiction with the observed cool core / non-cool core 
dichotomy. 



There are several outstanding issues which need to be 
addressed in more detail with follow up work. First we need 
to extend our study to both lower-mass galaxy groups and 
the most-massive superclusters to verify that our conjecture 
regarding the different impact of early AGN pre-heating on 
groups and clusters holds. At the same time, large num- 
ber statistics are mandatory to support our view that the 
two types of entropy profiles (cool and non-cool cores) are 
simply linked to different evolutionary stages of cluster as- 
sembly histories. It also remains to be shown that we can 
retrieve this result when metals are included in the simula- 
tions. Finally, the metallicities of both the simulated ICM 
and central massive galaxy under-estimate observations so 
far. It will therefore be important to quantify how resolu- 
tion, stellar winds and type la SNe alter metallicity profiles, 
cooling times, and, thus, the gas properties of the ICM. The 
effect of X-ray photo-ionisation of ion species which act as 
coolants of the proto-cluster intra-halo gas must also be ex- 
plored as it is possible that it will offset the large amount of 
extra metal cooling. In summary, we find that the cooling 
core structure is extremely sensitive both to the physics and 
the AGN feedback modelling, and consequently requires a 
more careful investigation. 
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APPENDIX A: VOLUME- WEIGHTED VERSUS 
X-RAY- WEIGHTED PROFILES 

A commonly adopted way to measure ICM gas properties 
is to weight them by their X-ray emission (oc n^VT). In- 
deed, as observations of ICM gas are often done in X-ray 
(although the Sunyaev-Zeldovich (SZ) effect ca n also be a 
good probe of the gas tempe rature, density, etc.. iKav et al.l 
l2008l : IProkhorov et all I2OI0I ). it seems natural for one to 
mass- weight or emission- weight one-dimensional gas profiles 
when comparing simulations with real data. However, the 
problem with such weightings is that strong density con- 
trasts will mask the contribution of diffuse gas. This typi- 
cally happens when a cold, gaseous galactic disc is present in 
the centre of a galaxy cluster. The central parts of the pro- 
files are dominated by the ISM emission because of its large 
density, and so the temperature will shift towards the peak of 
radiative cooling around 10^-10* K. One solution often used 
is to excise this ISM dominated emission by removing all gas 
elements colder than some temperature (~ 0.1 keV), which 
corresponds to the lower temperature limit of X-ray instru- 
ments {Rosat, Chandra, XMM-Newton). Another form of 
weighting which might provide a better co mparison to ob- 
serva tions is spectroscopic-like weighting (jMazzotta et al.l 
I2QO4I ). In this appendix, we discuss why none of these pre- 
scriptions is entirely satisfactory. 

In fig. lAll we compare spectroscopic-like weighting 
(first column). X-ray weighting (second column) and vol- 
ume weighting (third column) of entropy, gas density, and 
temperature profiles for the AGNOFFrun. The cluster in 
this simulation does not harbour any cold ISM component 
at z = 0.58 (recall it is restarted from the 2 — 0.58 out- 
put of AGNJETrun), and then very quickly develops a cold 
and massive galactic disc in the centre (see section 13.51 for 
details). The first thing to note is that the slopes and the 
levels of entropy, density, and temperature, are very insensi- 
tive to the weighting method sufficiently far away from the 
centre (r > 80 kpc), independent of redshift. Furthermore, 
when no cold gas disc is present (e.g. at z = 0.58, z — 0.50, 
and z = 0.42) there exist no significant difference between 
the different weighting methods all the way to the centre of 
the cluster. 

It is only when the cold ISM has built up significantly 
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Figure Al. A comparison of spectroscopic- like weigliting (first column), X-ray weigliting (second column) and volume weighting (third 
column) of entropy, gas density, and temperature profiles for the AGNOFFrun. Cells with temperature below T = 0.3 keV are excised 
in the profiles using spectroscopic-like weighting, while those with temperature below T = 1 keV arc excised for both the X-ray and 
volume weighted profiles. Colors correspond to profiles measured at different redshifts listed in the panel of the entropy profiles for the 
spectroscopic-like weighting case. 



that the profiles start to differ strongly. Even when cold 
gas cells (defined as cells with temperature T < 0.3 or 
T < 1 keV) are excluded from the analysis, X-ray and 
spectroscopic-like profiles are dominated by the disc com- 
ponent. The reason for this is that even though most of the 
cold gas in galaxies is at very low temperatures (T <^ 0.1 
keV), the gas at the interface between the disc and the ICM 
lies close to a very strong discontinuity, and as such is ex- 
tremely sensitive to any kind (numerical or physical) of dif- 
fusion process that mixes the two phases together, especially 
because it is extremely difficult to resolve the cooling length 
(Acooi = Cstcooi) . This disc 'skin' is extremely difficult to 
peel from the profiles except by applying a prohibitive den- 
sity/temperature cut-off. In addition to these 'skins', large 
galactic discs also sport a ffaring component at large radii, 
which is dense (p ~ 0.1 cm"''), hot (T ~ 1 keV ), and low 
entropy {K ^ 100 keV.cm'^). As X-ray emission is propor- 
tional to p^, this means that even when suppressing cold 
gas cells (T < 1 keV), profiles will still be dominated by the 
disc tail component. The characteristic signature of such a 
disk tail is apparent in the X-ray weighted profiles of fig. lAll 
(middle column) with cells below a temperature of 0.3 keV 
removed. At ~ 40 kpc from the centre, a dip is visible in both 
the entropy and temperature profiles, and a bump shows up 
in the density profiles. In contrast, volume-weighting the 
profiles, allows us to mostly probe the volume filling ICM 
without being strongly affected by the disc component of the 
central galaxy. For this reason, we adopted volume- weighted 
quantities throughout this paper. 
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